Self-Organized Bottleneck in Energy Relaxation 
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We study an energy relaxation process after many degrees of freedom are excited in a Hamiltonian 
system with a large number of degrees of freedom. Bottlenecks of relaxation, where relaxations of 
the excited elements are drastically slowed down, are discovered. By defining an internal state for 
the excited degrees of freedom, it is shown that the drastic slowing down occurs when the internal 
state is in a critical state. The relaxation dynamics brings the internal state into the critical state, 
and the critical bottleneck of relaxation is self-organized. Relevance of our result to relaxation 
phenomena in condensed matters or large molecules is briefly discussed. 
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Relaxation process to equilibrium has gathered much 
attention of physicists over several decades. Although the 
relaxation near equilibrium has been well formulated the- 
oretically, the process far from equilibrium is not yet fully 
understood, while relaxation phenomena after strong ex- 
citation can show novel and interesting behaviors. For 
example, single-molecule Myosin, a motor protein, shows 
anomalously long energy storage when biochemically ex- 
cited by ATP 1\. Such slow relaxations are also found 
in the dynamics of single enzymes 0, 1^ ■ 

When a system is excited weakly, it is described by 
a superposition of elementary excitations with weak in- 
teractions. When the excitation is strong, on the other 
hand, it is inevitable to take into account strong inter- 
actions of the excited modes, which may behave coop- 
eratively. Even though a general theory for such case 
may not be available, it is important to find a class of 
novel relaxation phenomena, and propose novel physical 
concepts associated with it. 

In this Letter, we report a relaxation process when 
quite a few degrees of freedom are highly excited, in a 
simple Hamiltonian system with large degrees of free- 
dom. The excited elements which interact strongly form 
a partial thermodynamic system, to which an internal 
state (with effective temperature etc.) is assigned. The 
relaxation depends on this internal state, and, in turn, 
the internal state dynamically changes with the relax- 
ation. We will show that, due to the interplay between 
the relaxation and the internal state dynamics, the par- 
tial system self-organizes a critical state, which forms a 
bottleneck in the relaxation course. We discuss its mech- 
anism, in particular by emphasizing distinction of it from 
critical slowing down and self-organized criticality j^. 

As a specific model we adopt a Hamiltonian system 
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The equations of motion are given by, 
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where 'N pendula are globally coupled, interacting each 
other through phase difference 27r(0i ~ ^i)- Each pendu- 
lum can have two types of motion; rotation at a higher 
energy and libration at a lower energy. 

This system is known to exhibit continuous phase 
transition in equilibrium in the thermodynamic limit 

^ oo Order parameter M for this transition is 



defined as, Me^"^""^ = YI^^^i e'^'""' 1^- When total en- 
ergy E is very small, M ~ 1, where all pendula li- 
brate. As the total energy increase, M decrease, where 
some pendula librate and others rotate, changing their 
types of motion temporally. Above the transition energy 
= 0.75/(27r)2, M ~ 1///V, and all pendula rotate 
almost freely. 

Energy relaxation when a single pendulum is excited 
is found to be very slow 0, with the relaxation time in- 
creasing as a stretched exponential form of the excited 
energy, as is discussed in terms of Boltzmann- Jeans con- 
jecture 0,0. 

Here, we are interested in relaxation properties when 
a part of the system composed of many pendula are ex- 
cited. At f = 0, pendula are simultaneously kicked, 
i.e. excited instantaneously, with the same momentum 
P^^^ . Here, we consider the case that the rest system 
is huge enough, and is regarded to be in thermal equi- 
librium. Instead of carrying out numerical simulations 
of the Hamiltonian system 0J with such huge degrees of 
freedom, we add a heat bath to the rest of the system; 
i.e. equations of the motion of the rest N —N'^^^ pendula 
are given by the Langevin equation. 
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FIG. 1: A typical time series of momenta of pendula. — 
64, T = 0.01,'iV(^' = 32, P'-'^ = 0.6. 
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FIG. 2: A typical time series of (a) the population of excited 
pendula A'^' , (b) the effective temperature of the excited part 
r'^' (©), and (c) the effective critical temperature T^^^' O- 
The data are the same as Figure Q 



where S,i(t) is white-gaussian random noise, that is 
= and (6(io)0(*o + 1)) - 2jTd,jS{t). 

A typical time series of the energy relaxation after the 
kick is drawn in Figure ^ The excited pendula that 
rotate relax one by one to librate with the rest of the 
pendula, so that the population of excited pendula 
decrease from to zero, whose time series is drawn in 
Figure 12- (a). It shows several plateaus, where the relax- 
ation is drastically slowed down. We call these plateaus 
as bottlenecks of the relaxation. 

With regards to this intermittent appearance of the re- 
laxation bottlenecks, we ask whether the bottleneck ap- 
pears randomly or follows some rule. If the latter is the 
case, the state of the excited partial system should be im- 
portant, since other huge degrees are kept at a constant 
temperature during the relaxation. 

Here, the interaction between the rotating and librat- 
ing pendula is small on the average, because the tempo- 
ral average of the term sm2Tr{9i{t) — 9j{t)) is negligible 
if the momenta of the two pendula are much different. 
Indeed, this is the origin of long relaxation time scale 
in the Boltzmann- Jeans conjecture 0. The interaction 



between two excited pendula, on the other hand, can be 
strong, if the momenta are not much different. In fact, by 
taking a Galilean frame moving with the average veloc- 
ity of the excited pendula, they often librate each other, 
keeping attracting interactions. 

Hence our system is represented by two partial sys- 
tems which are weakly coupled. One system is composed 
of excited pendula whose momenta ~ P^^^' , while the rest 
system has momenta ~ 0. Since the interaction between 
the two partial systems are small, it is relevant to con- 
sider a thermodynamic state of each partial system at 
the first-order approximation. The state of the relaxed 
part is basically determined by the heat bath with the 
temperature T. The state of the excited part, on the 
other hand, changes as the relaxation process progresses. 
Thus it is important to study a thermodynamic state of 
the excited part as an internal state of the whole system. 

As a quantity representing the internal state, we intro- 
duce effective temperature of the excited part as 

i 

by taking an incrtial frame with the center-of-mass mo- 
mentum of the excited part = Ylii^'^Pil ■^^^^ where 
Yj'^ denotes the summation over excited pendula 0|. 

To see how the relaxation depends on the internal 
state, we first define the first relaxation time r as the time 
when the first one of N'^^^ excited pendula loses energy to 
librate. In other words, r is the first escape time of a pen- 
dulum from the excited to the relaxed part. The relation 
between t and T*^^^ is plotted in Figure |31 where T*^^^ is 
time-averaged for a while just after the pendula are ex- 
cited. A peak is discernible. The temperature for this 
peak corresponds to the critical point for the phase tran- 
sition of the excited part tJ^^ = T^N^^^/N when 
the excited part is regarded as in equilibrium. 

We have also computed the time series of T^^-* and 
T^^^\ and plotted in Fi gureEKb) and (c), respectively. 
The plateau of N'--^^ actually corresponds to the time 
when the effective temperature of the excited pendula 
is around the critical temperature. From these results it 
is concluded that the relaxation bottleneck occurs when 
the internal state is around a critical point. 

One might think that this phenomenon is nothing but 
a critical slowing down of the excited system. This is not 
the case. First note that the whole system is not nec- 
essarily in a critical state. In fact we verified that the 
coincidence of this peak position to the effective critical 
point is qualitatively independent of the heat bath tem- 
perature T (not shown). Next note that such long relax- 
ation bottleneck is not found when disregarding the rest 
system. In fact we estimated the first escape time of a 
pendulum in only a single system, by calculating the time 
until one of the pendula first takes a momentum smaller 
than a threshold. This escape time increases with the 
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FIG. 3: The first relaxation time r versus the effective tem- 
perature of excited part T'^'. t'^' is computed from the time 
average of ^ over a period of 1000 just after the pendula are 
excited [ij, where the time scale of libration is the order of 
0.1 to 1. Instead of simple kick, we here attach two systems 
at t — 0, whose center-of-mass momenta are and 0. By 
varying the temperature at t < 0, T'^' is controlled. A'^ = 32, 
T = 0.01, iV<^' = 16, P'^ = 0.8. 2000 samples are computed 
in all. r is averaged for the samples within the bin size of 0.1 
for log scale. Tc'^' is the effective critical temperature defined 
as Tc*^) =TcN^'^'>/N O. 



decrease of the temperature, without showing any peak 
(not shown). Hence the present slow relaxation is inter- 
system relaxation, not mfra-system relaxation to which 
the so-called critical slowing down belongs. 

Then how is the excited system sticked to the criti- 
cal point through the interaction with the rest system? 
To answer this question, we consider how the internal 
state changes before and after one of the excited pen- 
dula relaxes. To see this change, we introduce effective 
temperature change AT^^ defined as, 



J\J{E) _ l'^ after -^before^ 



-T, 



(6) 



where T^lf^Q^g and Tj^^^^ is the effective temperature be- 
fore and after the first relaxation of an excited pendu- 
lum, respectively, while the factor N^^^ /{N^^ — 1) comes 
from the contribution of scaling The relation be- 

tween and Tjjfj^^g is plotted in Figure ^ where 

T^]ore ^'^'^ T^Tfter ^'^s time avcragcs over a given pe- 
riod just before and after r, respectively. The average 
of AT(^) is positive when T^^f]^^^ < T^^, while negative 



when T^fj^^g 



< This result indicates that the in- 

ternal state tends to be attracted into the critical state, 
through the interaction between the excited system and 
the rest. 

To analyze the above temperature dependence of 
AT^^, we note two factors for it; the increase of the 
temperature by work given to the excited part from the 
relaxed part, and the decrease of the temperature by an 
escape of a pendulum from the excited part. 



FIG. 4: 



versus Tj^gj^^^, computed from the simulations 



of Figure 01 T^^efore ^^"^ ^i/ter computed from the time 
average over a period of 1000 just before and after r [TH . re- 
spectively. AT^^' is averaged for the samples within the bin 
size of 0.1 for log scale. T^'^' is the effective critical temper- 
ature. 



First, a work is given from the relaxed part to the 
excited part, which causes temperature increase as the 
second law of thermodynamics indicates. Temperature 
change by the work is given by AT^^ = J dt T^^Kt) with 



f(E) 



(7) 



Using the order parameter of the partial systems 
^(*)g»0<*> ^ Y,/*^^'''^ /N'^*\ where * is £: or 7? [l^, we 
obtain the following expression, which is closed only by 
macroscopic variables: 



dt N 



dt V 
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Here E^^^ is the effective total energy of the excited part. 
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N 1 
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. The first term 



given by E^^^ — 2^ 

of the r.h.s. of (jSJ is a fast varying term, and can be 
neglected. With this approximation, the total energy 
change E^ of the excited part is proportional to M'^^ . 
Since AT^^^ ~ AE'^^/C^^, where C^^) is effective spe- 
cific heat of the excited part, we obtain. 



(9) 



Its temperature dependence in the thermodynamic limit 
of the excited part is drawn in Figure [S]- (a), which is 
positive at T^^ < tJ^ and vanishes for T^^) > T^'^^K 

Second, the temperature of the excited part decreases 
when a pendulum escapes from the excited part, which 
has the largest negative momentum at the center-of- 
mass inertial frame of the excited part. This factor 
is estimated by calculating the decrease of temperature 
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FIG. 5: Temperature dependences of {a)MC~^ and (b)— 
in the thermodynamic limit N ^ oo, which are obtained 
analyticaUy by considering canonical distribution In the 
present mean-field model, the critical exponent of M at T = 
Tc — is 1/2, and so the specific heat C does not diverge but 
takes a largest finite value, (c) is an example of combination 
of the two factor (a) and (b) with the ratio of 1:0.4. The 
temperature of horizontal axis is for a single system, and has 
to be scaled with N''^'^ /N to compare with Figure 0] 



expected to be general, as long as the system can exhibit 
continuous phase transition, and one part of the system 
is highly excited. 

The self-organized bottleneck at a critical state may be 
observed in condensed matter, and in molecules with suf- 
ficiently large degrees of freedom, when a partial system 
is highly excited. One candidate for the experimental 
verification is partial excitation of semi-conductor ma- 
terial by laser radiation, where so-called exciton Mott 
transition is found 13]. Also proteins show drastically 
slow relaxation when they are functioning as motors [l| 
or enzymes |3, |^ • The present bottleneck may shed new 
light on the understanding of such phenomena. 

The authors are grateful to S. Sasa and N. Nakagawa 
for discussion, and A. Shimizu for informing Ref |l,1| . 
This work was supported by a Grant-in- Aid for Scientific 
Research from the Ministry of Education, Science, and 
Culture of Japan. 



with the decrease of energy, by neglecting the interac- 
tion with the relaxed part. The escape of a pendulum, 
with momentum p^sc measured from the center-of-mass 
incrtial frame, takes away the energy A_Eesc = Ptsd'^ 
from the excited part. Since temperature is a function 
of the total energy as T^^) = T^^(E^^), the effective 
temperature change © can be described as, AT^^ — 
iV(^)/(Ar(s) _ l)T(E)(E(^)) - T(^){E(^) _ AE,,,). In the 
limit oo, this is estimated by the specific heat 

as, 

Ar(^) ~ -C7(^)"'A£;e,c oc ~C^^~\ (10) 

whose temperature dependence is drawn in Figure [3-(b). 

By combining above the two factors. Figure 01 can be 
explained qualitatively, as shown in Figure [Sl-(c). AT^^^ 
is positive at a temperature below Tc^^\ As the tem- 
perature approaches Tc^^' from below, AT^^ decreases, 
Above tJ-^\ there is only the negative contribution. 
Hence AT^^) chan ges its sign from positive to negative as 
the temperature is increased from below to above tJ^\ 

In summary, by using simple Hamiltonian system, we 
have discovered a bottleneck of energy relaxation, given 
by a critical state of an internal state of the excited part. 
This critical state is self-organized through the course 
of relaxation. The long-term relaxation from highly ex- 
cited state to equilibrium consists of alternation between 
approach to critical state and collapse of the state. In 
contrast to the so-called self-organized criticality the 
critical state is spontaneously formed without continu- 
ous driving, once a part of a system is highly excited. 
Although we have demonstrated this bottleneck critical 
state by a simple Hamiltonian system, the mechanism is 
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